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We study the Josephson current in a junction comprising two superconductors linked by a strong 
ferromagnet in presence of impurities. We focus on a regime where the electron (and hole) motion 
is ballistic over the exchange length and diffusive on the scale of the weak link length. The current- 
phase relation is obtained for both two- and three dimensional ferromagnetic weak links. In the clean 
limit, the possibility of temperature-induced 0-n transitions is demonstrated while the corresponding 
critical current versus temperature dependences are also studied. 



I. INTRODUCTION. 

The Josephson effect is a striking manifestation of 
quantum mechanics at macroscopic scales [1]. When 
a small current / is driven through a superconduc- 
tor/insulator/superconductor junction, no voltage drop 
occurs along the junction while a finite phase difference 
X appears between the two superconducting order param- 
eters of the leads. When the applied current exceeds a 
maximal (critical) value /c, a finite voltage appears across 
the barrier yielding a time-dependence of the phase. In 
many cases the stationary current phase relation (CPR) 
is well approximated by its first harmonic /(x) — h sinx 
and then the critical current is simply given by Ic — Ii- 
Nevertheless, theory gives room to higher harmonics, in 
particular at low temperatures. In fact, the only general 
requirement is that I{x) must be a 27r-periodic and odd 
function of the phase difference when time-reversal in- 
variance is respected. Thus the CPR may be expressed 
as a Fourier sum I{x) — '^„i ^rnsinmx where the coef- 
ficients Im are related to processes whereby m = 1,2, .. 
Cooper pairs are transferred through the weak link [2] . 

If the junction is inserted in a superconducting loop, 
the supercurrent is controlled by the applied magnetic 
flux <& which is directly related to the superconducting 
phase difference by x = 27r$/<I>o, 'i'o — h/2e being the su- 
perconducting flux quantum. In the absence of any bias- 
current or magnetic flux, the equilibrium state of a tun- 
nel Josephson junction usually corresponds to x = and 
/ = 0. In contrast, when the tunnel barrier contains mag- 
netic impurities, it was predicted that the phase differ- 
ence may be equal to tt at equilibrium. This may enable 
a spontaneous persistent current to flow in a loop com- 
prising a Josephson junction in such a 7r-state [3, 4]. This 
TT-shift is related to processes whereby electrons change 
their spin projection when passing through the insulating 
layer [5, 6]. Unfortunately this kind of 7r-state, generated 
by magnetic impurities in a host insulating layer, was 
never observed experimentally. In contrast it was further 
predicted [7, 8] and observed [9, 10] that a superconduc- 
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tor/ferromagnetic metal/superconductor (SFS) junction 
also exhibits transitions between zero and 7r-groundstates 
when the exchange energy h and/or the length L of the 
ferromagnet is varied [7]. The corresponding current- 
phase relation (CPR) /(x) were analysed for pure [7] and 
dirty [11] ferromagnets using respectively Eilenberger 
and Usadel equations [12-15]. In both cases, the criti- 
cal current of a SFS junction oscillates and decays when 
increasing the length or the exchange energy of the ferro- 
magnet. The oscillations of Ic{L) originate directly from 
the exchange interaction which induces a finite mismatch 
between the Fermi wavevectors of spin up and down elec- 
trons. Besides these oscillations, scattering by magnetic 
and nonmagnetic disorder strongly suppresses the criti- 
cal current when L is increased. In diffusive ferromag- 
nets, this overall decay is exponential on the typical scale 
^1 = y^D/h, D being the diffusion constant. In contrast, 
in the pure limit, a finite Josephson current may be ob- 
served up to much larger length scales on the order of 
= '^f/T, vp and T being respectively the Fermi veloc- 
ity in the ferromagnet and the temperature. In particu- 
lar at zero temperature, the decay becomes a power law, 
namely Ic ~ in the case of a three dimensional pure 
ferromagnetic weak link [7]. In the absence of disorder, 
this critical current suppression results from the super- 
position of many distinct single channel CPRs associated 
with independent transverse channels. Accordingly this 
decay is expected to be less severe in low dimensional fer- 
romagnets for it corresponds to angular averaging over 
quasiclassical trajectories with different angles with re- 
spect to the junction axis. 

The first evidence of the 7r-state, in SFS junctions, 
came with the observation of the nonmonotonic behavior 
of the critical current as a function of temperature [9] and 
of the ferromagnet length [10]. The weak ferromagnetic 
alloys chosen for these pioneer experiments enable to ob- 
serve the transitions in relatively large junctions within 
the ten nanometers scale. Further support was provided 
by magnetic diffraction patterns of DC squids comprising 
a SFS TT-junction in one arm and a usual tunnel junction 
on the other arm [16]. Furthermore spontaneous persis- 
tent currents were reported in a loop interrupted by a tt 
junction [17] and imaged in arrays of tt- Josephson junc- 
tions [18]. Finally, multiple transitions between zero and 
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FIG. 1: SFS junction. The contact area are far smaller 
than the typical superconducting lead (S'l,S'_r) transverse 
size. The metallic ferromagnet (F) in the middle part can 
be either two- or three- dimensional. 



TT-groundstates were observed by varying the ferromag- 
netic layer thickness of a SFS junction [19, 20]. 

Experimentally obtaining the CPR is much more dif- 
ficult than simply measuring the critical current. Only 
recently a few CPR experiments were implemented suc- 
cessfully in the case of SFS junctions [21] and SNS junc- 
tions [22-24]. In fact at sufficiently low temperature, a 
small second harmonic of the CPR is always present both 
in SNS and SFS junctions, but it is usually completely 
eclipsed by the large amplitude of the first harmonic. The 
SFS Josephson junctions are a natural playground to ob- 
serve unambigously the second harmonic since the first 
one vanishes at the zero-7r transition. In particular, the 
second harmonic was detected as a tiny minimum super- 
current at the crossover between the zero and tt states, 
and also revealed by the related Shapiro steps [25] . In the 
highly transparent limit, the issue of the sign and mag- 
nitude of the second harmonic was adressed in presence 
of uniaxial [26] and isotropic [27] magnetic scattering in 
the ferromagnet. Moreover finite transparency or weak 
interfacial disorder may also modify the second-harmonic 
[28]. 

Zero-TT transitions were also observed in smaller junc- 
tions comprising strong ferromagnets like Fe, Co, Ni or 
permalloy [20, 29, 30]. These novel experiments are 
performed in an interesting regime which differs both 
from the pure clean or dirty limits extensively studied 
so far. Owing to the extremcUy large exchange energy, 
the period of the critical current oscillations ^2 = vp/h 
is smaller than the mean free path £ = vpT while the 
ferromagnetic bridge is still longer than the mean free 
path: ^2 ^ £ ^ L. Bergeret et al. investigated the 
first harmonic of the Josephson current in this particular 
regime hr ^ 1 [31]. A theoretical analysis of this second 
harmonic in this particular regime is still lacking while it 
was already detected experimentally [20]. 

Up to now the physics of the 7r-state and second har- 
monic were mostly investigated in the three-dimensional 
case due to the lack of lower-dimensional ferromagnets. 
Recenlty novel systems like graphene or thin films of mag- 
netic semiconductors became available as promising can- 
didates for realizing two-dimensional SFS junctions. For 



instance coating graphene with Pd may produce itinerant 
magnetism [32] while alcaline coating is likely to produce 
superconductivity [33] . Hence tailoring a graphene sheet 
with appropriate metals on top may induce SFS hetero- 
j unctions within the carbon atoms plane. At the present 
time, only two experiments have been reported on SNS 
junctions made with graphene [34, 35], while O-tt tran- 
sitions were predicted in graphene based SFS junctions 
[36]. Another experiment which has triggered the inter- 
est on two-dimensional SFS junctions is the measurement 
of a supercurrent through a long, L = 0.3 — l/xm, half- 
metallic ferromagnet chromium oxide (Cr02 ) film [37] 
wherein singlet superconductivity should be destroyed on 
a much smaller length scale. The spin singlet to triplet 
conversion (at the interfaces) was proposed to explain 
the finite Josephson current. It was shown that triplet 
correlations penetrate a ferromagnet on much longer dis- 
tances [38] . Another possibility is that the surface of the 
film is less spin-polarized than the bulk or even antifer- 
romagnetically ordered. In this latter scenario, a singlet 
supercurrent may bypass the half-metallic ferromagnet 
by flowing within a two-dimensional SFS surface junc- 
tion. 

In this paper, we study both two- and three- 
dimensional SFS junctions with strong exchange field and 
moderate disorder, namely in the limit hr ^ 1. Us- 
ing Eilcnberger equations and perturbative expansion in 
l/hr, we obtain the current phase relation and in partic- 
ular its second-harmonic at the O-tt transition which was 
recently observed in the three-dimensional case [20]. In 
the pure limit, we show that the two-dimensional criti- 
cal current is suppressed as L~^^^ instead of in the 
three-dimensional case. Accordingly we suggest the pos- 
sibility to observe enhanced critical current in planar SFS 
junctions made of magnetic semiconductors or graphene 
with an induced ferromagnetic order. 

After introducing the formalism in Sec. II, we in- 
vestigate the pure limit in Sec. Ill with special em- 
phasis on the two-dimensional case and on temperature- 
induced O-TT transitions. In Sec IV, the CPR of two- and 
three-dimensional SFS junctions are obtained in the limit 
/ir 1 in relation with experiments [20]. 



II. SFS MODEL AND FORMALISM. 

We study a superconductor-ferromagnetic- 
superconductor (SFS) Josephson junction in the 
geometry represented in Fig.l. The ferromagnet consists 

in a single ferromagnetic domain characterized by its 
exchange energy h, length L = Lx and transverse 
dimension(s) Ly (and L^)- We assume that the super- 
conducting order parameter is Aoe*-*^/^ (rcsp. Aoe^''^/^) 
in the left (right) superconducting lead while the Fermi 
velocity vp is the same everywhere. The contacts 
between superconductors and the ferromagnet are 
completely transparent and spin inactive. Besides the 
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geometrical parameters of the junction, three typical 
lengths are of primary importance for the Josephson 
effect. On the one hand the superconducting coherence 
length ^0 = Uf/ttAq and the ferromagnetic exchange 
length = vp/h are related to the strengh of the 
superconducting and ferromagnetic order parameters 
respectively. On the other hand disorder is characterized 
by the elastic mean free path I = vpT where r is 
the average time between impurity scattering events. 
Henceforth, we adopt units with h = Ub = ^■ 

Previous theoretical studies were mostly performed in 
the diffusive limit I ^ ^i?, L or in the pure limit £ ^ (_p, L 
using respectively the Usadel and the Eilenberger equa- 
tions (without self-energy terms due to disorder) [12]. 
Recent experiments have opened to possibility to inves- 
tigate the regime < ^ < i [20, 29, 30]. This later 
regime cannot be described by the Usadel equation since 
the disorder induced self-energy is no longer the domi- 
nant energy. It is thus necessary to use the Eilenberger 
formalism [31] and the exchange energy as the large pa- 
rameter which enable perturbative expansions. 

In our simple model, the quasi-classical Eilenberger 
Green functions g = g^{x,Vx), f = fu{x,Vx) and 
= {x,Vx) depend only on the center-of-mass co- 
ordinate X along the junction axis Ox, and on the angle 
of the quasiclassical trajectories with respect to Ox. 

In the ferromagnetic weak link, < L/2, the Eilen- 
berger equations read 

f vAg={2T)-\f{f+)-f+ (/)), 

vAf = -2{uj + ih) + -{g{f)-f (g)) , 
^ 1 

-vM+ = -2 (c^ + th) f+ + -{g (/+) - /+ {g)) . 

T 

(1) 

Here lo = ttT (2n + 1) are the Matsubara frequencies and 
Vx — vp cos 9 is the Fermi velocity vector [12]. The brack- 
ets denote averaging over the Fermi surface. 

In the superconducting leads, which are assumed to be 
clean, the Eilenberger equations read : 



vAg = A*/ - A/+ 
Vxdxf = -2a;/ + 2Ag 
vAf+ = 2wf+ - 2A*g 



(2) 



with A = Aoe^^/"^ (resp. Aqc-'^^/^) for the left (resp. 
right) electrode respectively. In the whole paper. Aq = 

is the temperature de- 
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pendent superconducting gap. 

In the limit /ir ^ 1 studied thorough this paper, one 
may assumes as a starting approximation that (/) = 
(/+) = 0. Then demanding the continuity of the gen- 
eral solutions of Eqs. (1,2) at the interfaces x = ±L/2 
yields the quasiclassical Green functions over the whole 
junction. In particular, in the ferromagnet, |a;| < L/2, 
one finds that the normal Green function 



g{x,Vx) 



A2 
^0 



sinh$ 



is independent of the position x. The upper (lower) sign 
of the denominator corresponds to the positive (negative) 
sign in the velocity projection. We have defined = 
+ Aq and the effective phase 
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(4) 



which contains all the relevant parameters of the junc- 
tion. In experiments, the exchange field h is always larger 

than the temperature, and thus the contribution ujL/vx 
may be safely neglected in the above expression for 

The supercurrcnt density is given by the following qua- 
siclassical expression [12] 



j (x) = 2TTev'^^^T ^ {v, Img^ {v,)) 



(5) 



where the temperature T is expressed in energy units and 
v^'^ is the density of state at the Fermi level per spin and 
per unit volume (resp. surface) for d = 3 (resp. d = 2). 
The corresponding current / (x) is obtained as the flux 
of i (x) through a section of the weak link. 

Finally the groundstate energy of the junction can 
be deduced by integrating the CPR according to the gen- 
eral formula 



2tt dE^ 



(6) 



LO 



Q. Q, uj sinh $ ± cosh $ 



(3) 



where $0 = h/2e is the superconducting flux quantum. 
The — TT phase transition occurs when the X = and 
X = TT groundstates are degenerate, namely when Eq = 

E^. 



III. SFS JUNCTION IN THE PURE LIMIT. 

In this section, we consider the pure limit L ^ £ with 
special emphasis on the two-dimensional case. Indeed, 
the three-dimensional and the one-dimensional cases are 
well known for both small [7] and large [39, 40] exchange 
energies. After briefly recalling the single channel results, 
we obtain that the low temperature critical current of 
a two-dimensional Josephson junction decays as L~^l'^ , 
namely more slowly than the dependence charac- 
terizing three-dimensional ballistic weak links. We ob- 
tain the CPR and study the second harmonic at the O-tt 
transition, where the first harmonic cancels. We also 
study the possibility of O-tt transitions induced by vary- 
ing the temperature at a given length. Finally, the Ic^) 
curves are obtained and compared to recent experiments 
[19, 20]. 

A. Single channel case. 

State of the art ferromagnetic wires typically still con- 
tain a large number of transverse channels. Nevertheless, 



4 



the analysis of the single channel case is both necessary 
and instructive since, in the ballistic limit, it is the build- 
ing block for evaluating the multichannel supercurrent in 
higher dimensions. When considering a single transverse 
channel, the angular averaging (..) in Eq.(5) reduces to 
a discrete sum over 9 = and 9 = w which yields the 
following current-phase relation: 



o-=±l ^ 



— - cos — sm — , 

2T 2 2 



(7) 

where a = 2hL/vp and llf' = ev'^^'^ irvp^o- Using 
Eq.(6), one obtains the energy of the junction 



E, 



(1) 



<T=±1 



In 



, , Ao X + era 
cosh — - cos — - — 
2T 2 



(8) 



as a function of the phase difference %• The typical en- 



ergy scale is given by Eq^'' = 2T7ny^''' Hvp. The x = and 
TT groundstates are degenerate for regularly spaced values 
of the accumulated phase a, namely at a„ = 7r/2 -|- titt. 
These critical vahics of the parameter a can be reached 
by varying either the length L or the exchange energy 
h, whereas they are insensitive to temperature variation. 
Hence in an hypothetical single channel SFS junction, it 
would be impossible to drive the O-tt transition by vary- 
ing the temperature only. 

Close to the critical temperature, T « Tj,, the lin- 
earization of Eq.(7) yields a nearly sinusoidal current- 
phase relation / (x) = Ic{a) sinx whose critical current, 
Ic{a) = /Q^''(Ao/2rc) cosa, cancels for a = ■7r/2 -|- titt, 
namely at the O-tt transitions. Moreover the oscillatory 
behavior of Ida) is not damped when a, or equivalently 
L and/or h, are increased. This behaviour is due to the 
absence of angular averaging in the single channel situa- 
tion. 

At lower temperatures, T <C Tc, the CPR becomes 
nonsinusoidal. The corresponding critical current Ic (a) 
is obtained numerically by maximazing the current den- 
sity I (x) given by Eq.(7). This critical current Ida) 
also exhibits periodic oscillations as a function of a = 
2hL/vp. In contrast to the situation for T Ri Tc, the cur- 
rent Ic (a) is finite at the cusps owing to the presence of a 
sizeable second harmonic. Finally, it is very instructive to 
check the sign of this second harmonic I2 {an ) at the zero- 
TT transitions where first harmonic cancels /i(q!„) = 0, 
since this sign is related to the order of the transition. 
It turns out that hian) > which yields a discontinu- 
ous (first-order) phase transition between the to the tt 
phases. Otherwise, namely for l2{an) < 0, the transi- 
tion would have been continuous (second-order) with a 
groundstate corresponding to an arbitrary value of the 
phase difference, distinct from or tt [14]. 



B. Two-dimensional case. 

In a two-dimensional SFS junction, the supercurrent is 
the sum of the currents carried by independent transverse 
modes. Using the angular averaging (..) = J d9/2'K{..) 
appropriate for planar junctions, one obtains the follow- 
ing CPR: 



a^dy 



X tanh — - cos sm (9 

V2r 2 J 2 ^ ' 

(2) (2) 

where a = 2hL/vF and /q = ev^ ttvf^o- The corre- 
sponding critical current Ic (a) is shown in Fig. 2 for sev- 
eral temperatures. In contrast to the single channel case, 
the oscillations of Ida) are damped due to the angular 
averaging over many transverse channels having each a 
distinct CPR. For T <Tc, the curves Ida) exhibit cusps 
where the critical current remains finite instead of the 
cancellations observed for T « Tc. 

The zero- and 7r-groundstate energies Eo{a) and -Ex (a) 
depend both of a and temperature according to: 



Eo{a) 



E, 



(3) 
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E, 



(3) 
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cosh ( ^ cos - 
2T 2 



cosh I ^ sm - 



a^dy 



(10) 



a'^dy 



(11) 

The values of the parameter a = 2hL/vp where the 
(O-tt) transitions occur are obtained by solving Eoia) = 
E^(a) which may be rewritten as 
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, , Ao y 

cosh COS — 

2T 2 



dy 



cosh 



2T 2 



(12) 



At low temperature T < Tc, we have checked that the 
solutions adT) for Eq.(12) coincide with the location of 
the cusps in the Ic (a) curves. Moreover, a.c{T) decreases 
when the temperature is lowered, sec Fig. 2. The phase 
diagram in the a-T plane is similar than the three dimen- 
sional phase diagram shown in Fig. 4. It should be also 
emphasized that this transition is not accompanied by a 
non-monotonic behavior of the Ic {T) curves in contrast 
to the dirty case [19]. 

Near the critical temperature Tc, the current-phase re- 
lation is nearly sinusoidal with a first harmonic given by 



hia)=I^ 



_ r(2) Ao 



dy- 



a cos y 



"^TcJa y^^y^ - a^ 

and a second harmonic given by: 

o? cos 2y 



dy 



(13) 



(14) 
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FIG. 2: Critical current as a function of a = 2hL/vF in a 
two-dimensional SFS junction. The temperature is increased 
from the upper curve (T « 0) to the lower one {T xTc). The 
— TT phase transitions occur at the kinks of the Ic(<y) curve. 
For a given — tt transition, the positions of a„{T) of these 
kinks exhibit weak temperature dependence. 



FIG. 3: Critical current Ic (a) versus a = 2hL/vp for T ^ Tc 
, from Eq.(13) in the two-dimensional case (solid curve) and 
from Eq.(20) in the three-dimensional case (dashed curve). 
In the two-dimensional case, the period of Ic is larger than 
in the three-dimensional case. The overall decay of Ic is also 
much slower in the two dimensional case. 



Both /i(a) and Ii{a) exhibit an oscillatory dependence 
on a, see Fig. 3. In the limit a ^ 1, the asympotic be- 
haviors are given respectively by 



Ii {a, T Tc) Ag \/tt cos a — sin a 



r(2) 



2Tr_ 2 



/a 



and by: 



7-(2) 

-'o 



Aq \ ^ \/7r cos 2a — sin 2a 



2Tc 



24 



2a 



(15) 



(16) 



Thus in the regime a 3> 1, the O-tt transitions occur for 
the values a = 7r/4 + nn where the first harmonic of the 
CPR cancels. Moreover the second harmonic is positive 
at these transitions. 



C. Three-dimensional CEise. 

Finally we consider the well-known CPR for three- 
dimensional SFS junctions [7] : 



a^dy 

y3 



X tanh — - cos sm (17) 

2T 2 J 2 ^ ' 

in order to compare with the two-dimensional results 
described in the previous paragraph. Here /q ' = 

(3) (3) (3) 

bvq LyLzTTVp^o and Eq = 2X^1^^ Hvp- The zero- and 



TT— groundstate energies Eo{a) and E^{a) depend both 
of a and temperature according to: 



Epja) 
p(3) 

E^{a) 



E, 



(3) 





In 



In 



cosh 



Ao y 

-T^ COS - 

2T 2 



, , Ao . y 
cosh — - sm - 
2T 2 



a^dy 

yd 

a^dy 

y3 



(18) 
(19) 



By solving numerically Eq (a) = Eq (a) at arbitrary tem- 
perature, we obtain the curves ac{T) shown in Fig.4. In 
principle, the system can experience — tt phase transi- 
tion by lowering the temperature at some given a, pro- 
vided this value of a is close to a critical value. It is 
important to note that the range wherein the (0 — tt) 
phase transition can take place by tuning the tempera- 
ture is smaller for the three-dimensional case than in the 
two-dimensional one. 

Near the critical temperature Tg, the first harmonic is 
given by 



T I \ r(3) Ao /""^ 
h(a) =^0 J cosy 



a^dy 

y3 



(20) 



which is plotted on Fig. 3. The asympotic behaviors are 
given by 



h{a,T^n)^-I, ——- 
in the limit a ^ l.The second harmonic 
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FIG. 4: Phase diagram in the a-T plane showing the 
transition hnes between the and the tt phases for a 
three dimensional SFS junction. We have used Ao (T) = 

1, 764Tc tanh (l, 74y/T/Tc - Ij . 

also exhibits an oscillatory dependence with respect to 
a. Thus in the regime a 3> 1 and close to Tc, the O-tt 
transitions occur for the values Q!„ = mr where the first 
harmonic of the CPR cancels. The second harmonic is 
positive at those points, I2{n7r) > 0, indicating first order 
transitions. 



D. Temperature dependence Ic{T) 

For appropriate values of the ferromagnetic layer 
length and exchange energy (corresponding to a = a„), 
it is possible to pass through the 0— tt transition point by 
changing the temperature, as shown in Figs. 4 and 5 for 
a = 3.72. This kind of temperature induced — tt tran- 
sition was actually achieved in experiments using dirty 
weakly ferromagnetic alloys [19]. It was observed that the 
critical current exhibits a nonmonotonic T-depcndence 
with a cancellation at the — tt transition. Moreover the 
TT state can be either the low (e.g. at the first node) or 
the high temperature phase (at second node) [19]. In con- 
trast, only monotonic variations of Ic(T) were reported 
in experiments with strong ferromagnets in the clean (or 
moderately dirty) limit [20, 29]. 

Here we have obtained that, in the pure limit, the crit- 
ical current does decrease monotonously when the tem- 
perature is increased, as shown in Fig. 5. Nevertheless 
the temperature dependence of the critical current I,, (T) 
exhibits very distinct shapes (e.g. in Fig. 5 for a = 3 and 
a = 3.72) being either convexe, concave or almost lin- 
ear depending on the value of a = 2Lh/vF- Although 
wc have shown extreme cases (in Fig. 5 for a = 3 and 
a = 3.72), most of our IciT) curves are almost linear 



FIG. 5: Critical current Ic (T) versus temperature for a three- 
dimensional junction. Four different values of a = 2hL/vF are 
shown: a = 2.25 (solid line), a = 3.0 (dotted line), a = 3.72 
(solid line) and a = 4, 5 (dashed line). The critical current 
decreases monotonically when the temperature is increased. 
The dot on the a = 3.72 curve indicates the (O-tt) transi- 
tion, which is not associated witli any singular beliavior in 
the Ic (r) curve. In the inset, the above values of a are indi- 
cated on the Ic (a, T Tc) curve. 



in agreement with the experimental curves reported in 
[20, 29]. For future experiments, we suggest the observa- 
tion of the concave curves (e.g. a = 3.72 in Fig. 5) as a 
signature of the O-tt transition. The corresponding mea- 
sure is quite challenging since it corresponds to the tem- 
perature dependence of a minimum of the critical current 
which is given only by the contributions of higher har- 
monics (to > 2). Similar features were predicted in the 
limit of large exchange fields using Bogoliubov-de Gennes 
formalism [40]. Here we confirm that this change in the 
concavity of the Ic{T) curve near a O-tt transition is still 
predicted in the limit of moderate exchange fields in com- 
parison to the Fermi energy. In contrast, the Ic{T) curves 
are non monotonic in the dirty limit when the junction 
passes the O-tt transition. We explain this discrepancy 
by the fact that the critical current at the O-tt transi- 
tion vanishes in the dirty limit whereas it is still finite in 
the pure limit (due to the important contribution of high 
harmonics). 



IV. SFS JUNCTION WITH IMPURITIES. 

Experiments on SFS junctions comprising dilute mag- 
netic alloys are correctly described within the Usadel 
equation framework [19, 26-28, 41], because the exchange 
energy is smaller than the disorder level broadening, 
namely r/i <C 1, and far smaller than the Fermi energy. 
In this regime the electron (and hole) motion is diffu- 
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sive with a mean free path smaller than both and 
L. Recently experiments were performed in the oppo- 
site regime, rh > 1, using strong ferromagnets, like Fe, 
Co, Ni or Permalloy. Then the electron (and hole) mo- 
tion is ballistic over the ferromagnetic length scale ^p, 
while being still diffusive on the scale of the weak link 
length L. In particular this situation implies that the 
parameter a — 2L/^p is very large. The first harmonic 
of the CPR for three dimensional weak links has been al- 
ready found [8, 31] by solving the Eilenberger equations 
for large rh. In this section, we calculate the amplitudes 
of the first and second harmonic both in two- and three- 
dimensional SFS junctions. The analytical expressions 
obtained here can be used as a starting point to inter- 
pret the three-dimensional experiments by Robinson et 
al. [20] and future investigations on graphene-based SFS 
junctions [36]. 



-0.01 
-0 . 02 
-0 . 03 




FIG. 6: First (red curve) and second (blue curve) harmonic 
as functions of L/l for a three-dimensional junction. We have 
chosen 1/£,f ~ 5. The corresponding analytical expressions 
Eqs. (23,24) are derived in the regime ^ I '!S L. Hence the 
lower limit for L/l is one. 



A. Three dimensional case. 

Here we investigate the CPR I{x) — h sin x + 
/2sin2x -t- ... of a three dimensional junction. We have 
obtained the first harmonic 



A3) 



uj>0 



and the second harmonic 



A3 

^0 



(w + ny 



L 



(23) 



Rc{jEiiz)-E3{2z) 



(24) 

where Iq^^ — ei^Q^^ LyLzTrvpAa, LyLz being the junction 
area, and z — L/l + ia ~ L/l + 2iL/£^p. The functions 
Ei{z) are defined in the appendix. The first harmonic 
corresponds to the result previously obtained in [8, 31]. 
The second harmonic is usually smaller than the first 
one, except when Ii cancels. Then we obtain that the 
magnitude of I2 is perfectly measurable though small. 
Moreover, the sign of I2 when /i = is very instructive 
since it determines the order of the transition. When the 
first harmonic cancels, the second harmonic is finite and 
positive. This is the usual case where one observes a dis- 
continuous jump of the junction between the zero and tt 
states. For negative I2 when Li = 0, one would expect a 
continuous transition and the realisation of a, ip junction 
[14]. For the experimentally relevant regime of moder- 
ate L/l, we always obtain that the second harmonic is 
positive at the transition (see Fig. 6). For larger val- 
ues of L/l, we have also observed negative I2 although 
numerical artefact cannot be excluded. 

We now study how the supercurrent is suppressed 
when the weak link length or the exchange field is in- 
creased. The asymptotic behaviors of the harmonics at 
a ^ 1 are given by 



-757= . Re 



ui>0 



{UJ + 17)2 



L/1 + 2iL/ij. 



(25) 



and 



h_ 

r(3) 

■'0 uj>0 



^ ^ 8Age-2V' ( (L/l - 2iL/ip)e-^^^/i^ 



ir i'^ + ^y V {L/l + 22L/^p)^ 



(26) 

Besides the exponential suppression Li ^ e (and 
Q-2L/1 j^-j^ ^Yxe real parts in Eq.(25,26) provide 
damped oscillations as a function of a, as shown in Fig. 
7. 

The sum over Matsubara frequencies yields the tem- 
perature dependence of the harmonics. At low tempera- 
ture T — *■ 0, the sum over Matsubara frequencies becomes 
an integral that can be done analytically yielding: 



A3) 



■Re 



L/l + 2iL/^F 



-L/l 



and 



/2 _ 8 ^ fiL/l-2iL/^p)e-^^^/i^ 



,-2L/l 



(27) 



, (28) 



= — Re , 

j(3) 15 V iL/l + 2iL/^p)^ 

at large a. 



B. Two dimensional case. 



For planar junctions, both the first and second har- 
monics were unknown in the strong ferromagnet regime 
£,P I L. Using the procedure described in appendix, 
we have evaluted these harmonics as 



i6 4'^7^E^^^i^^^2(z), 



(U + 0)2 



(29) 
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FIG. 7: First (red curve) and second (blue curve) harmonic 
as functions of L/Z for a two-dimensional junction. We have 
chosen 1/£,f = 5. The corresponding analytical expressions 
Eqs. (29,30) are derived in the regime ^ I ^ L. Hence the 
lower limit for L/Z is one. 



and 



requires weak exchange fields in practice. The overall 
decay or damping of the critical current as a function 
of the length/exchange field of the ferromagnet is much 
slower (L~^/^) in the two-dimensional case than in the 
three dimensional case {L~^). Moreover, the shape of the 
critical current versus temperature curves changes when 
one closely approaches a O-tt transition. Hence for fu- 
ture experiments, we suggest the challenging measure of 
the temperature dependence of a minimum of the crit- 
ical current, which is essentially given by the contribu- 
tions of higher harmonics (m ^ 2). The corresponding 
/c(T) should differ markedly form the usual almost linear 
curves obtained so far away from those minima [20, 29]. 

We have obtained the current phase relation for SFS 
junctions comprising strong ferromagnets in the presence 
of moderate disorder, namely in the limit rh > 1. We 
have calculated the second harmonic, in particular at the 
O-TT transition, for both two- and three- dimensional SFS 
junctions. In the three dimensional case, we have com- 
pared our result with recent experiments performed in 
the regime rh > 1 [20]. 



A3 



Re 



F^iz)-F,{2z: 



Here /, 



(2) 



,(2) 



- {to + nr 

(30) 

TTVp^oLy, Ly being the width of the 
planar weak link, and z — L/l + 2iL/S^p. The functions 
Fi{z) are defined in the appendix. From Fig. 7, one ob- 
seves that the second harmonic is finite and positive when 
the first harmonic cancels. 



VI. APPENDIX. 



In this appendix, we derive the current-phase relation 
in the regime L 3> ^ ^ ^f, or equivalently a rh 1. 
We start with the normal quasiclassical Green function 
which is uniform in the ferromagnet: 



C. Comparison with experiments 

In the experiments [20] , the values of the parameter hr 
are respectively hr = 3 for NigoFe2o, 2.8 for Co, 1.62 for 
Fe and 0.5 for Ni, using the parameters (/i, vp, I) provided 
in [20]. The experiments [29] performed on Nb-Ni-Nb 
junctions correspond to /ir ~ 1. Hence these experiments 
cover the onset of the regime hr > 1, or equivalently < 
I. The oscillations of the critical current are reported for 
weak link lengths L not exceeding few mean free paths 
I, othewise the signal would be too small (due to the 
exponential suppression by the factor e~^/'). The period 
of the oscillations is approximatively 0.5 in units of L/l 
in agreement with our results for hr — 3, see Fig. 7. 
Quantitative comparison between our theory and these 
experiments are hindered by the fact that band structure 
effects and interface quality may strongly influence the 
magnitude of the Josephson current. 



V. CONCLUSION. 

In the absence of impurities, we have demonstated that 
temperature-induced O-tt transitions are possible both 
in two and three dimensional SFS junctions though it 



^0 



sinh$(6l) 



Q Lu smh<P{e)±n cosh <i>{e)' 



where +/— corresponds to the sign of cos 9 and 



2m = ix 



{9)L/l- 



cost 



(31) 



(32) 



with a — 2hL/vp. We first consider the case of positive 
Matsubara frequency w > 0. Then, one expands gu, {6) 
in powers of X\ — exp [— 2A<I>(6')] as 



9^ (^) = 1 - 



2A§ 



2A4 



{uj + 



xi 



(33) 



where A is the sign of cos 6. The modulus of X\ is a small 
parameter because \X\\ ~ exp{—2X{g)L/{lcos6)) with 
(g) > (In first approximation (g) = sgn{ijj)) and L ^ I. 
This expansion defines a self-consistent problem since X\ 
contains {g) which in turn is evaluated using X\. As a 
first iteration, we evaluate g^^ (9) using (g) = 1 in the 
expression of Xx. Then we perform the angular averaging 
to obtain the first order correction to {g) . Denoting {g) = 
1 + .91 ( Iffil ^ 1) O'^G obtains, in the three dimensional 
case: 
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where: 



91 = - 



2Agcosx 
{u + 0)2 

2Agcosx 



7r/2 

/ rf6'sin6»exp ( ^) (34) 

J V cost'/ 



i^2(2) 



CSC 



J. 

^35^ The current is related to the average: 



(41) 



with z = L/l + ia and En{z) = J dy y "e 

1 

The current is related to the quantity: 

7r/2 
A 2 f 

{v.g^ m = -J^^ff^ J d9 sine cose{X^ - Xq) 



(36) 



{cu + n)4 



7r/2 

J desmecose{xl- xl) 



In this expression, the first integral (proportional to Ag) 
must be evaluated within the approximation: 



- = -2isinxexp ( ^-- \ (l-fi'i-^) 

V cosC/ \ cosy J 

whereas the zero-order approximation {g) = 1 : 

/ 2z 

X'^ — Xq = —2i sin 2x exp 



cos^ 



(37) 
(38) 



is sufBcicnt for the second integral (proportional to Aq). 
Substituing Eq. (37,38) in Eq.(36) yields: 



{Vx9u, {0)) 



(w + 17)2 



E3{z) sinx 



+ 



2ivFAt ^i.£2(^)_^^(2.))sin2x 



(w + 0)4 V / 



and finally Eq. (23,24). 

In the two dimensional case, the first order correction 
to {g) reads: 

7r/2 

2A^cosx /■ / L/l + ia\ 



7r(a; + 0)= 
2Ag cos X 

" (W + 0)2 



COS^ 



(40) 



it/2 

{v.gu. {6)) =^F cose[g{e) - gin - 6)] 



tt/2 



+ 



2vfAI fde 

^-^y-cos^(Xe-Xe) 


Tr/2 

, / -COS«-X|). 

(W + 0)4 y TT ^ ® ^' 



1x12 

2vfA% f de 



Proceeding to the same approximations as for the three- 
dimensional case, one obtains: 



, 4isinxi'_FAn ^ , , 



+ 



(W + 0)2 

4isin2xuFAo / L 



[UJ + 0)4 \l 



-Fiiz)-F,{2z) 



and finally Eqs. (29,30). 
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